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Abstract. We study the linear post-Newtonian approximation to general rela- 
tivity known as gravitoelectromagnetism (GEM); in particular, we examine the 
similarities and differences between GEM and electrodynamics. Notwithstanding 
some significant differences between them, we find that a special nonstationary 
metric in GEM can be employed to show explicitly that it is possible to intro- 
duce gravitational induction within GEM in close analogy with Faraday's law of 
induction and Lenz's law in electrodynamics. Some of the physical implications 
of gravitational induction are briefly discussed. 



PACS numbers: 04.20.Cv 

1. Introduction 

In electromagnetism, the combined dynamics of charged particles and electromagnetic 
field are consistently described by Maxwell's field equations and the Lorentz force law. 
However, in the linear perturbation approach to gravitoelectromagnetism (GEM), one 
recovers the Maxwell equations for the GEM field, but the corresponding Lorentz force 
is recovered, to first order in v/c, only when we deal with a stationary GEM field. 
This explains why some authors (see, for instance, [TJ [5]) have treated GEM only 
for stationary fields and the issue of existence of gravitational induction in analogy 
with Faraday's law of induction is therefore absent in such treatments; moreover, it 
has been argued recently that in general relativity such an analogy does not even 
exist 3 . On the other hand, time- varying GEM fields have been implicitly considered 
by many authors (see, for instance, 0] — HD])- In fact, some gravitational Faraday 
experiments were proposed in [4] based on the existence of gravitational induction in 
analogy with electrodynamics. The purpose of the present paper is to show explicitly 
that general relativity does indeed contain induction effects; these turn out to be, 
despite the differences that have been mentioned, on the whole closely analogous to 
electromagnetic induction effects. 

In general, GEM covers those aspects of general relativity that can be best 
explained via an electromagnetic analogy. In this paper, we work mainly within 
the linear GEM scheme; therefore, it is necessary to review briefly the relevant 
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aspects of this linear post-Newtonian approximation to general relativity [B] in order 
to render the present paper essentially self-contained. Consider the curved spacetime 
generated by a localized slowly rotating "nonrelativistic" astronomical source. In 
the linear approximation, the spacetime metric can be written as = t]^ + h^, 
where ?7 M „ is the Minkowski metric tensor with signature +2 in our convention and 
h^v is a first-order perturbation. Under a slight transformation of the background 
coordinates x M = (ct, x), x M i— > x M — e' 1 , the gravitational potentials h^ v transform 
as h,j, v i— > h^y + e^ v + e u>fl . Henceforth, the potentials are considered to be gauge 
dependent, while the background global inertial coordinate system is in effect fixed. 
The spacetime curvature is, however, gauge invariant. It is useful to introduce the 
trace-reversed potentials = h fiv —-zhri lill with h = tr(/i M „). Imposing the transverse 
gauge condition h^" \ v = 0, the gravitational field equations take the form 

- 16vrG , . 

□v = — — T »» ■ 

The general solution of (jl.ip is given by the special retarded solution 

h»- C 4 y | x _ x /| d * . i 1 -^ 

plus a general solution of the homogeneous wave equation that we simply ignore in 
this work. In the linear GEM approach, all terms of <3(c -4 ) are neglected in the metric 
tensor. It then follows from Eq. (|1.2| that for the sources under consideration here 
/ioo = 4<&/c 2 , hoi = —2Ai/c 2 and hij = 0(c~ 4 ), where $(t,x) is the gravitoelectric 
potential and A(t, x) is the gravitomagnetic vector potential. The spacetime metric 
is thus given by 

ds 2 = -c 2 ^1 - 2^J dt 2 - ^(A • dx)dt + (\ + 2^J ^diW , (1.3) 

where far from the source the dominant contributions to the GEM potentials can be 
expressed as 

, GM . GJxx 

$ = , A= — . (1.4) 

r c r'> 

Here M and J are the inertial mass and angular momentum of the source, r = |x|, 
r ^> GM/c 2 and r ^S> J/(Mc). Let us note that the gauge condition implies that 

^$ + V-Qa^=0. (1.5) 

This is related to the conservation of mass-energy of the source via Eq. (|1.1[) . That is, 
let T 00 = pc 2 and T° l = cj l , where j M = (cp, j) is the mass-energy current of the source; 
then, Eq. (|1.5p is equivalent to j 7 '. A , = 0. It is possible to define the gravitoelectric 
field E and the gravitomagnetic field B in close analogy with electrodynamics 

E = -V3> - -^dt Qa^ , B = V x A . (1.6) 
It follows from these definitions that 

VxE=-i 9t (V), V.(V)=0, (1.7) 
while the gravitational field equations (jl.ip imply 

V-E = 47rGp, V x [ ^B ) = i<9 t E + — j. (1.8) 

V 2 / c c 
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These are the Maxwell equations for the GEM field. The particular form of these 
equations is based on a special convention [TT] that makes it possible to employ 
the standard results of classical electrodynamics in the GEM framework. This is 
accomplished by assuming that the source has gravitoelectric charge Qe = GM 
and gravitomagnetic charge Qb = 2GM. Moreover, a test particle of inertial mass 
to has gravitoelectric charge qE = —m and gravitomagnetic charge qs = —2m in 
this convention. The signs of (qE,qB) are opposite to those of (Qe,Qb) due to the 
attractive nature of gravity; furthermore, the ratio of gravitomagnetic charge to the 
gravitoelectric charge is always 2, as the linear approximation of general relativity 
involves a spin-2 field. This circumstance is consistent with the fact that the ratio of 
the magnetic charge to the electric charge of a particle is unity in Maxwell's spin-1 
electrodynamics. We note that the magnetic charge employed here is different from 
the magnetic monopole strength, which is always strictly zero throughout this work. 

Given Maxwell's equations for the electromagnetic field, Faraday's law of 
induction simply follows, for instance, from the consideration of the temporal variation 
of the magnetic flux linking a static closed circuit. A similar approach in the GEM 
case would fail, however, as the line integral of E along the closed circuit does not in 
general correspond to the work done by the gravitational field of the source. This is 
the crucial point and to clarify the situation, it is therefore necessary to investigate 
the motion of a free test particle in the linear GEM scheme. 

We must now discuss the analogue of the Lorentz force law in our linear GEM 
framework. The geodesic equation for the motion of a free test particle is 

du^ 

~~dT 

where t/c is the proper time and w M = dx^/dr is the unit four- velocity vector of the 
test particle. The Christoffel symbols are given by 

c 2 r° 0Al = - $ „, c 2 r% = 2A (lJ) + o , (i.io) 

c 2 r 00 = - $ - 2A 4 , , c 2 r l Q3 = o + e t]k B k , (1.11) 
c 2 T % ok = fe + - 6 jk $ ti . (1.12) 

The geodesic equation can be reduced via it^ = 7(1, f$) with v/c to 

r(h (l-^ 2 )$, + 2^[$ )i -A (i>j) ^], (1.13) 



T% a uP U ° = , (1.9) 



7 dt 
dv 1 



(l+(3 2 )$ ti -2(0xB) i +2A i 



- /3*(3 - /? 2 )$ o + 2/?/^[% fe) /3 fc - 2$,,] . (1.14) 
Moreover, u^u^ = —1 implies that 



- = l-P 2 --(l + P 2 )<f + -0-A. (1.15) 



- 2 =l-[3 2 -\{l + (i 2 )<S> + - 
7 Z c c 

For a stationary source (<9t$ = and 9 t A = 0), Eq. (|1.14p reduces to 

dv v 

m— — = — mE — 2m— x B, (1.16) 
dt c y ' 

when velocity-dependent terms of order higher that (3 — v/c are neglected. In the 

case of a general nonstationary source, however, the equation of motion (|1.14j) does 

not correspond to the Lorentz force law and this implies that the electromotive force 

does not in general have a simple analogue in GEM. 
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Though the gravitational analogue of the Lorentz force law has a more 
complicated form in GEM, we intend to show via a special nonstationary GEM 
metric that induction effects can still exist in close analogy with electrodynamics. 
The motivation for our approach comes from a detailed consideration of the 
gravitomagnetic clock effect. This is briefly discussed in the next section. 



2. A nonstationary GEM metric 

We start our analysis with a brief discussion of the gravitomagnetic clock effect, since 
there is an important heuristic connection between gravitational induction and this 
effect. Consider circular equatorial geodesies about a Kerr source of mass M and 
angular momentum J. The exterior spacetime metric is given by 

ds 2 = - c 2 di 2 + - (dp 2 + X dl9 2 ) + (p 2 + a 2 ) sin 2 

X 

(cdt- asin 2 6»d0) 2 , (2.1) 



2Mp f l( 2 



where M = GM/c 2 , a = J/(Mc) > is the specific angular momentum of the Kerr 
source and 

£ = p 2 + a 2 cos 2 9, x = P 2 ~ 2Mp + a 2 , (2.2) 

in Boycr-Lindquist coordinates. The geodesic equation for a circular equatorial orbit 
reduces to 

At I a , . 

aq> u>k c 

where lok is the Keplerian frequency, ujk — (GMj p 3 ) 1 / 2 , for the orbit with fixed 
"radius" p > 2M and 9 = n/2. The upper (lower) sign in Eq. (|2.3|) refers to a 
co-rotating (counter-rotating) orbit with respect to the sense of rotation of the Kerr 
source. It follows from (12.311 that 



27r a , 
t± = ± 2tt- , (2.4) 

ujk c 

where t + (t_) is the period of prograde (retrograde) circular motion in terms of the 
proper time of the static inertial observers that are infinitely far from the source. We 
are interested in physical situations where the orbital motion is far from the source, 
i.e. p ^> 2M and p » J/(cM), so that geodesic motion is possible in opposite 
directions for the same orbital "radius" p. Then t+ — £_ = 4irJ/(Mc 2 ) illustrates the 
gravitomagnetic clock effect. This remarkable result, which is independent of G and p, 
holds to lowest order for the proper times of clocks in orbit around the source as well. 
It has already been discussed in a number of papers (see, for example, [51 ITS"]); 
therefore, we concentrate here on the fact that the prograde motion is slower than 
the retrograde motion. Specifically, let v± — 2irp/t± be the relevant speed of motion 
according to the static inertial observers at spatial infinity (p — > oo); then, to first 
order in au^fc <C 1, 

GJ 

v±~v K T -5-0 > ( 2 - 5 ) 
c z p z 

where Vk = P^k- 

Imagine an ensemble of identical Kerr spacetimes except for different magnitudes 
of J. As J increases in this ensemble, v+ decreases and V- increases. Let us first 
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note that this circumstance cannot be interpreted as "inertial induction" [14] . as the 
effect is simply the opposite of what such a Machian interpretation would predict - 
we return to this subject in section 5. On the other hand, if one could turn this 
kinematic situation into a dynamic one in terms of the temporal variation of J, then 
one could at least heuristically interpret the change in speeds in terms of induced 
currents due to a time-varying flux of the gravitomagnetic held. To this end, we need 
a solution of the held equations of general relativity that would correspond to a Kerr 
solution but with a time- varying J. It turns out that such a solution exists [15] . but 
only in a rather approximate form within the linear GEM framework. We therefore 
turn to a description of this solution (see also Appendix A) . 

Let us hrst note that the Kerr solution can be put into the form of metric (| 1 . 3|) 
with potentials ([1 .4|) once the isotropic radial coordinate r 

(2.6) 

is introduced in Eq. (12. ip and the resulting metric is linearized in M/r and a/r with 
x = r sin 9 cos <f>, y — r sin 9 sin cj> and z — r cos 9. 

Consider next a spacetime metric of the GEM form (| 1 . 3[) with potentials 

GM G , .Jxx , nr . 

*= , A=- J + Jit)— =- , 2.7 

r c r 6 

where the magnitude of the proper angular momentum of the source varies linearly 
with time, i.e. J(t) = (Jo + Jit)J. Henceforth we assume that J = z and Jo > 0. 
Substituting potentials (|2.7p in the GEM held equation (jl.ip and gauge condition 
(jl.5p . we hnd that the latter is satisfied and the former gives the effective source of 
the spacetime metric in the form 

Too = Mc 2 5(x), r 0l = ^c[J(t) x V],;<5(x) . (2.8) 
It follows from the dynamical equations for the source 

= (2.9) 

that 



T — - — 

v ~ 8^ 



(j x V)iVj- + ( J x V)-,Vi- 



(2.10) 



where J = dJ/di = JiJ is a constant vector. While the mass-energy current (|2.8p 
is confined to the origin of spatial coordinates, the stresses (|2.10p are distributed 
throughout space and fall off as r~ 3 for r — » oo; however, this unusual circumstance 
has no impact on the viability of this nonstationary GEM spacetime, since T{j is 
independent of time and Eq. (|1.2p then implies that hij = 0(c~ 4 ) for the stresses 
(|2.10p . Thus this special source satisfies the requirements of the linear GEM scheme 
and hence generates an acceptable GEM field. That is, metric (|1.3p together with 
potentials (|2.7p represents a solution of the linearized field equations that is valid at 
the first post-Newtonian order of approximation. Further details about this time- 
varying but nonradiative solution are given in [15] and Appendix A. It is clear that 
the linear perturbation approach eventually breaks down over a sufficiently extended 
period of time due to the linear temporal variation of the gravitomagnetic potential in 
Eq. (|2.7p . On the other hand, we wish to avoid any complications associated with the 
instants of time at which the temporal variation is switched on and off. We therefore 
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consider a certain time interval after the temporal variation is switched on and before 
it is switched off such that our linear GEM approach is valid; in fact, we will always 
work within this interval of time for which 2|A| <§; c 2 . 

There exist radiative solutions of Einstein's field equations of the type originally 
due to Vaidya in which the mass (and hence possibly angular momentum) of the source 
can vary with time due to the emission or absorption of radiation. Such solutions are 
not of interest here. Instead, we concentrate on nonradiative solutions in which the 
angular momentum slowly varies with time. For instance, the Earth slowly loses 
angular momentum with time due mainly to tidal friction. Conservation of angular 
momentum implies that the orbital angular momentum of the Moon about the Earth 
increases in this process. The Earth-Moon distance thus slowly increases; the rate of 
orbital expansion at present is about 4 cm per year. The special nonstationary solution 
can be employed to discuss the physical implications of the temporal variability of the 
gravitomagnetic field. In the case of the GP-B experiment performed in orbit about 
the Earth, for example, the influence of such variability turns out to be negligible 
[15] . However, our considerations may be of interest in the astrophysics of rotating 
gravitationally collapsed configurations that exhibit variability. 

The next section is devoted to the illustration of gravitational induction and 
Lenz's law using the special nonstationary GEM solution (|2.7p . 

3. GEM induction and Lenz's law 

The special nonstationary spacetime given by potentials (|2.7p involves a static 
gravitoelectric field and a linearly time- varying gravitomagnetic field 

t 2c r 6 

B = — (J + Jit)i[3(J-x)x- J] . (3.2) 

c r J 

Thus the gravitational analogue of the displacement current is zero in this case. To 
first order in v/c, the analogue of the Lorentz force law is given by 

m ^ = _ m £_^ vxB , (3.3) 
at c 

where £ can be expressed as 

. GMx 2GJxx 

£ = — T^T~ ■ ( 3 - 4 ) 

The distinction between E and £ in general relativity is the root of the difference 
between the electromagnetic and gravitational inductions. A free particle initially at 
rest picks up an azimuthal speed due to the force term (2m/c)dtA in Eq. (|3.3|) . The 
general motion of a free test particle in this nonstationary spacetime is studied in 
some detail in the next section; moreover, certain aspects of the motion of spinning 
particles and light rays have been briefly considered in [T5] , 

Consider now a closed circuit in the equatorial plane of the variable rotating 
source as depicted in Figure 1. The circuit is bounded by an inner circle of radius r\ 
and an outer circle of radius r 2 . The radial parts of the circuit are infinitesimally close 
to each other. For the sake of concreteness, let the circuit consist of a perfect fluid 
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Figure 1. Schematic diagram of the closed equatorial annular loop for a thought 
experiment that illustrates gravitational induction. The loop is assumed to be 
sufficiently far from a nonstationary rotating source. 



that is at rest and fills a very thin tube. Once the circuit is placed in the time- varying 
gravitomagnetic field as in Figure 1, the flux of this field through the circuit is 

r =/B.dS = -^ (i )(I-l). (3.5) 

Moreover, we note that 

E . d< =I^(±-±) , (3.6) 

where d£ is an element of the circuit in the direction depicted in Figure f . Thus it 
follows from Eqs. ([53)1 and ([53)1 that 

E ' 1 '—hw' < 37 » 

which, as expected, is in accordance with Maxwell's equations for the GEM field. It 
turns out, however, that Eq. (|3 . 7|) is not the analogue of Faraday's law of induction 
in the gravitational case. 

In electrodynamics, the quantity evaluated in Eq. (|3.6p would be the 
electromotive force (emf), which is the amount of work done per unit electric charge. 
The gravitational analogue of this concept would be the amount of work done per 
unit gravitoelectric charge. We therefore define — on the basis of Eq. (|3 . 3[) — the 
gravitomotive force (gmf ) to be 

g= je-di, (3.8) 

where £ is given by Eq. (13. 4|) . since the v x B term does not contribute to the work. 
Thus the gravitomotive force is the circulation of £ rather than E. Calculating the 
gmf for the loop in Figure 1, we find 

^'I-IV (3.9) 

n r 2 J 
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Thus the gravitational analogue of Faraday's law of induction turns out to be 

"-If- 

Indeed, the actual induced current is four time larger than that given by Eq. (|3.7p , 
which is evident from the factor of four difference between the coefficients of 9<A 
in Eqs. p.4j) and (|3.ip . Nevertheless, these considerations show theoretically that 
gravitational induction exists and — within the linear GEM framework — is closely 
analogous to the Faraday law of induction in electrodynamics. 

It remains to show that the direction of the induced current is such as to 
oppose the change that caused it. Let us first remark that for a long straight line 
supporting a mass current /, Eq. (|1.8[) implies that the gravitomagnetic field has 
closed circular field lines of radius r around the line current as in electrodynamics; 
moreover, the magnitude of the field is 4GI/(cr) and its direction follows from the 
usual right-hand rule. In fact, this is the gravitational analogue of the Biot-Savart 
law of electrodynamics. For a moving test particle of mass m, the gravitoelectric 
charge is — m; hence, the direction of its current is opposite to its velocity. As the 
gravitomagnetic field increases (J > 0), the fluid (in the tube in Figure 1) begins to 
co-rotate from rest due to the induced azimuthal acceleration (2/c)dtA. The speed of 
this co-rotation is larger in the inner circle than in the outer circle of the loop, so that 
a clockwise motion develops in the loop corresponding to an induced counter-clockwise 
current and hence positive flux through the loop that opposes the increasingly negative 
flux of the gravitomagnetic field of the source. This is the gravitational analogue of 
Lenz's law of electrodynamics. 

Our treatment of GEM induction and Lenz's law has been based upon our specific 
thought experiment and the use of the special nonstationary solution. The main results 
are, however, quite general. Assuming <9 t $ = and working to first order in v/c, we 
find that the equation of motion (|3.3[) is generally valid with £ = — V$— (2/c)8j A and 
B = V x A for a nonstationary A(t, x). Thus with the definitions of the gravitomotive 
force Q and the gravitomagnetic flux J 7 , Eq. (13 . 1 0[) is generally satisfied. This is then 
the general form of the law of gravitational induction in GEM and the minus sign on 
the right-hand side of Eq. (|3.10|) is in accordance with Lenz's law. 

As demonstrated in section 2, in an ensemble of stationary Kerr spacetimes with 
varying J, the speeds of prograde (retrograde) circular equatorial orbits decrease 
(increase) with increasing J. It is interesting to demonstrate this effect dynamically 
using the nonstationary linearized Kerr spacetime. This is the purpose of the next 
section. 

4. Motion in a time-varying gravitomagnetic field 

We start with the equation of motion (|1.14jl adapted to potentials (|2.7|) . namely, 

dv GMx GM , , , „ , 2G j x x 2 „ 
— + — ^ = — [4(x • v)v - i/x + —— v x B 

dr r° err c r° c 

-^p[J.(xxv)](x.v)v (4.1) 

and inquire whether equatorial circular orbits are possible in this case. With x = 
r cos 0, y — r sin </> and z = 0, Eq. (|4.ip reduces to 

^ = - ( 4 - 2 ) 
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2 GM / v 2 \ 2GJ(t) 

where u = ro<^> and <p = d<f>/dt. Solving the quadratic equation for v to linear order in 
Newton's constant, differentiating the outcome with respect to time and comparing 
the result with Eq. (|4.2p . we find that Eqs. (|4.2|) - (|4.3[) are inconsistent so long as 
Ji 0. Thus no circular geodesic orbits exist in the equatorial plane of the source 
for Ji 7^ 0. On the other hand, for Ji = 0, we find from Eqs. (|4 2 [) -(|4.3 |) that for a 
circular orbit 



_ GM GJ 
VK± \l — ^ (44) 

in agreement with Eq. (|2.5p . Here the upper (lower) sign refers to a prograde 
(retrograde) orbit and terms beyond the linear order in Newton's gravitational 
constant have been neglected. Once J starts to increase linearly with time, the 
prograde (retrograde) equatorial circular orbits tend to spiral outward (inward); 
moreover, the azimuthal speeds of prograde (retrograde) orbits decrease (increase) 
according to Eq. (|4.1|) , as would be expected from our discussion of the 
gravitomagnetic clock effect in section 2. These results follow from the perturbative 
studies of Eq. (|4.1|) that are presented in the rest of this section. 
It is useful to write Eq. (14. 1|) as 

dv GMx ^ , A „, 

d^ + — =F ' (45) 
where F is the linear GEM relativistic perturbing acceleration. Let us digress here 
and note that to 0(c~ 2 ), F should also contain the nonlinear post-Newtonian term 
4G 2 Af 2 x/(c 2 r 4 ), which we have neglected in our linear treatment. For F = 0, the 
test particle follows a Keplerian orbit; in this case, the Newtonian energy En, orbital 
angular momentum Ljv and the Runge-Lenz vector Rjv (all per unit mass of the test 
particle) are conserved. These are given by 

1 , GM , s 
E N = , Ljv = x x v, 4.6 

2 r 

Rjr = v x Ljv - —x , (4.7) 

r 

so that in the presence of the perturbation F, we have 
dEN dLjv _ 



dt dt 
d~R N 



(4.8) 



= F x (x x v) + v x (x x F) . (4.9) 
dt 

For our present purpose, it suffices to concentrate on the energy equation and 
compute F • v using Eq. (|4.1|) along any unperturbed Keplerian circular orbit around 
the source. The unperturbed orbital plane can have an arbitrary inclination angle i 
as in Figure 2. For a circular orbit x- v = and hence the Newtonian energy equation 
reduces to 

d (\ 2 GM\ 2GL N ■ . 

- -i) 2 = „ , Jcosi , (4.10) 

dt \2 r J c 2 r 3 y 1 

where Ln — |Ljv|. Thus for Jcosi > 0, the Newtonian energy increases and hence 
the circular orbit tends to spiral outward, but for J cos? < 0, the orbit spirals inward; 
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moreover, for a polar orbit (cosz = 0), it remains unchanged within the orbital plane 
regardless of J. 

More generally, let us assume that the unperturbed orbit is an ellipse with 
semimajor axis Ro and eccentricity e in the {X, Y) plane. The correspondence between 
the (x,y,z) and (X,Y,Z) coordinate systems is illustrated in Figure 2. The ellipse 
can be represented by 

i?o(l-e 2 ) 



1 + e cos rj 



(4.11) 



where usq > is the corresponding Keplerian frequency, i.e. uJq = GM/Rq. 

Working in the {X, Y, Z) coordinate system, we introduce cylindrical coordinates 
(R, ip, Z), so that Eq. (|4.5p can be expressed in general as 

, GMR 

R-RV 2 + {R2 + Z2)3/2 =F R , (4.13) 
GMZ 

z + WTWr^ Fz > (4 - 15) 

where Fr and F v are given by 

Fr = F x cos ip + F Y sin ip, F v — —F x sin ip + F Y cos ip . (4.16) 

To compute the linear perturbation away from the Keplerian orbit due to 
relativistic effects, we evaluate Fr, F v and Fz along the unperturbed orbit and seek 
a solution of Eqs. (|4T3)) - (|4~T5)) of the form 

R = f(l + U), ip = r 1 + W 1 Z = rH, (4.17) 

where U, W and H are only considered to linear order. To simplify matters, it is 
convenient to express Eqs. ()4.13j) - (|4.15p in terms of the new independent variable 
rj instead of t, where £(77) is given by Eq. (|4.12[) . A lengthy, but straightforward 
calculation results in the linear perturbation equations [16] 



drj 2 drj 1 + e cos 77 

2U) = B , (4.19) 



d / dW 
dr/ \ dr/ 
d 2 H 



drj 2 



H = C , (4.20) 



where 



r 3 



(A,BX) = TT (F R ,F V ,F Z ) (4.21) 



J N 



and L 2 N = GMR (1 - e 2 ). We impose boundary conditions on Eqs. (|4. 18[) (14.20)) such 



J N 

that 



U = W = H = at 77 = %, (4.22) 

dU dW dH . . ... 
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Thus Eqs. l|3~Tg|t - P~Tg|) now take the form 
d 2 U l + 4ecos?7 



dr] 2 1 + e cos rj 

dW P 

^- + 2U = B( V ')d v > 



U = A(v) + 2 / B(n')dr/ , 



(4.24) 



(4.25) 



We must now evaluate A, B and C in order to solve the linear perturbation 
equations. The contribution of various source terms simply superpose; therefore, we 
limit our attention to the dominant relativistic terms up to linear order in v/c. Hence 



2GJxx 2 

-s ; v X B 



where B is given by Eq. (|3.2[) . It follows that with this F, we have 
2Gcosi J{t) 



A = 



B 



C = 



c 2 Lm f 
2Gcosi 



c l, n 
2G sini 

C L N 



_ ■ uj R esmri 

—rcosr/j + sm ^^ -f 3ecosn)J(t) 
V 1 — e 2 



(4.26) 

(4.27) 
(4.28) 
(4.29) 



where J = Ji, J(t) = Jo + Jit and 

J(t) = J + — (n - 2esin?7) + 0(e 2 ). (4.30) 

OJQ 

To simplify the solution of Eqs. (|4.18[1(|4.20[) . the source terms can be written as 
expansions in powers of the eccentricity as well as 

U = U + eUi + 0(e 2 ) (4.31) 

and similarly for W and H |16j . The boundary conditions would then apply term by 
term in such expansions. It can be shown that 
2 cos i 

A= ,, n { Jqlo + J\r\ + e[J uj cosr] + Ji(ricos?? — 2sinri)] + Q(e 2 )} , (4.32) 



B 



C 



Mc 2 
2 cosi 

Mc 2 
2 sini 

Mc 2 



{ Ji — e[J o uj sin 77 + Ji(?/sinr; + cos?/)] + 0(e 2 )} , 
{2 Jqujq sin i] + Ji (2rj sin 77 — cos 77) + e[3 Jqloq sin 77 cos 77 



(4.33) 



+Ji(3rysin7jcos?7 + cos 2 ry — 6sin 2 77)] + 0(e 2 )} . (4.34) 

Let us first consider the perturbation on an initially circular orbit. The 
expressions for Uo, Wo and Hq contain cumulative (secular) terms as well as harmonic 
terms. For instance, Uq is given by 
2 cos % 

Uq = [2Ji(T7-7jo) + Jom + Jit]-(J ujo + Jii]a) cos(r7-?7o)-3Ji sin(77 - 77o)] .(4.35) 

The dominant secular terms are given by 



6 cosi 



5 cosi 



-Jir? 2 



sin z 



rJi?7 cos?7 



(4.36) 



Mc 2 " u Mc 2 ' ' 1 Mc 2 ' 
Thus for Jx cos i > 0, the orbit spirals outward and the azimuthal velocity, given 
generally by 



dt ~ VT - ? 



dW 

(l + ecos??) ( l + E/~+ — ) > 



(4.37) 
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tends to decrease, since 

dW 4 cos i 

U ° + ^~-JMT J ^- (4-38) 

Next, we note that to first order in eccentricity the orbital perturbation can be 
calculated using Eqs. (|4.32[) - (|4.34p . The dominant secular terms turn out to be 

3cosz o . ... 6cosi 2 rj _ sini . 

1 ~ ~~MJ JlV Sm7? ' 1 ~~1Zc2" 1?? COS?7 ' 1 ~ ~M? Jir?Sm2r? • (4 - 39) 
In principle, one can continue this procedure and determine ?7, W and i? to all orders 
in eccentricity. 

In summary, the motion of a free test mass in the variable gravitomagnetic field of 
a central source is such that if the test particle starts from rest, it tends to move in the 
same sense as the source for J > and in the opposite sense for J < 0. On the other 
hand, if the test mass is already in almost periodic motion about the source, then 
for J > 0, the prograde (retrograde) motion tends to slow down (speed up) and the 
opposite takes place for J < 0. Our thought experiments illustrating the analogues of 
Faraday's law and Lenz's law in section 3 and appendix B involve test masses starting 
from rest. It has not been possible to provide simple thought experiments to illustrate 
in a similar way the behavior of test particles that are already in orbit about the 
central source. 



5. Discussion 

The purpose of this work has been to provide an explicit treatment of gravitational 
induction. The main ingredients of our discussion include the acceleration term 
(2/c)9fcA in the GEM force law — see for instance Eq. 1)1.140 — an d our special 
ansatz 1)2.70 for a gravitomagnetic vector potential A that varies linearly with time. 
It has been shown that, despite the existence of certain differences in the force law 
between GEM and electrodynamics, it is nevertheless possible to establish a certain 
analogy between gravitational induction and electromagnetic induction. 

The acceleration term (2/c)dtA in the gravitational force law has been 
traditionally interpreted to be in essence responsible for a Machian inductive action of 
accelerated masses such that a test mass would accelerate in the same direction as the 
acceleration of neighboring masses (see pages 100-103 of [14] ). However, our analysis 
of geodesic motion in the special nonstationary spacetime in this paper demonstrates 
explicitly that this Machian interpretation cannot be generally maintained in general 
relativity |17j . Indeed, for J > 0, nearly circular prograde orbits experience azimuthal 
deceleration rather than acceleration. 

Astronomical sources generally have variable angular momenta. For instance, 
external electromagnetic breaking torques tend to slow the rotation rates of pulsars 
with constant moments of inertia. The implications of our preliminary results for such 
systems are beyond the scope of this work. 
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Appendix A. Kerr metric with a = (ct 

The special nonstationary solution of linearized Einstein's equations has a nonzero 
Einstein tensor that modulo its symmetry can be expressed in spherical polar 
coordinates as 

G „ . (A.1) 

crr A 

away from the origin of spatial coordinates. Here J = J(t)z, J — J± and all the other 
components of the Einstein tensor vanish for r > 0. 

It is interesting to calculate the Einstein tensor for a Kerr metric obtained by 
changing a constant J to J{t) = Jo + Jit — or equivalently changing the constant 
parameter a to (Jo + J\t)/{Mc) — with Ji ^ in Eq. (|2.1[) . The result is quite 
complicated; to simplify matters, we make a time translation t t— ► t — Jo/Ji in the 
new time-dependent Kerr metric in Boyer-Lindquist coordinates. The result is metric 
(|2.1|) with a replaced by 

a = C,ct, (A.2) 

where £ is a constant dimensionless parameter given by ( — Ji/(Mc 2 ). If £ = 0, 
we recover the Schwarzschild metric for which the Einstein tensor vanishes. Thus an 
expansion of the Einstein tensor of this Kerr metric in powers of the small parameter 
£ should indicate how close this solution is to the exterior vacuum field of a stationary 
source. In this way, we find that Eq. (|2.1[) with a = £ci has an Einstein tensor such 
that modulo its symmetry G r g = 0, G tt = 0(( 4 ), while G tt j, and Gg^ are 0(( 3 ). 
Moreover, G tr , G t g, G rr , Ggg and G^ are all 0(C 2 ), but 

GM sin 2 6 , . „ , 

G ^ = - 3 c^r-2GM/c^ + °^- (A ' 3) 
To linear order in Newton's gravitational constant, Eq. (|A.3j) reduces to Eq. (jA.lj) . 



Appendix B. Induced current in a circular loop 

Imagine a circular loop of radius r consisting of a perfect fluid that completely fills a 
very thin tube and is at rest around a source of mass M and angular momentum J. 
The plane of the loop makes an angle of i with respect to the equatorial plane of the 
source as in Figure 2. 

We are interested in the current that is induced in the loop when J is linearly 
dependent upon time. The circular loop can be considered to be the boundary of a 
hemisphere of radius r above the plane of the loop; we use the surface of the hemisphere 
for the calculation of the flux of the gravitomagnetic field of the source through the 
loop. Using the spherical polar coordinates associated with the (X, Y, Z) system and 

z = cos i Z + sini Y , (B-1) 

the flux through the loop turns out to be 

T= —GJ{t) cos i . (B.2) 

cr 

Similarly, the gmf can be simply calculated from Eq. (|3.8[) . and we find 

0— £GJ«.i. (B.3) 
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X 

Figure Bl. Schematic plot of an unperturbed Keplerian orbit in the (X,Y) 
plane. The spatial part of the background global inertial frame is given by the 
(x, y, z) coordinate system. With respect to this, the (X, Y, Z) system is rotated 
by a constant inclination angle i about the x axis. 



Thus Eq. (f3~TU)) is satisfied by the results given in Eqs. (fR2|) and (|R3|) . 

It is clear from Eq. (|3.3[) that the fluid particles start from rest and move in the 
prograde sense for J > due to the presence of the acceleration term (2/c)<9tA. The 
corresponding induced current will be flowing in the retrograde sense and is such that 
its negative gravitomagnetic flux through the hemisphere opposes the increasingly 
positive flux of the source. The induced current is proportional to cosi, so that it is 
maximum for the loop in the equatorial plane and vanishes for a polar loop. 

It is important to remark that it is not possible to use the plane of the loop that 
passes through the source for the calculation of the flux. To see this, imagine instead of 
the surface of the hemisphere of radius r, the "sombrero" surface given by the annular 
region of inner radius r$ and outer radius r together with an upper hemisphere of 
radius rg that just avoids the source. The flux through each part of this surface is J- ' a 
and J 7 si respectively, where 

T A = — I - - — ] GJ(t) cos i, T s = —GJ(t) cos i . (B.4) 
c \r r s J cr s 

Note that for the effective source of the nonstationary solution, Eq. (|2.8|) . rs — > 0. 
Nevertheless, the net flux is T — Ta + J~s given by Eq. (|B.2|) . 
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